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Abstract 

A search has been made for lepton-flavor-violating interactions of the type 
e + p —>■ £X, where i denotes a fi or r with high transverse momentum, at a 
center-of-mass energy, y/s, of 300 GeV with an integrated luminosity of 47.7 pb _1 
using the ZEUS detector at HERA. No evidence was found for lepton-flavor vi- 
olation and constraints were derived on leptoquarks (LQs) that could mediate 
such interactions. For LQ masses below y/s, limits are set on X eq yJJkq-, where 
X eqi is the coupling of the LQ to an electron and a first-generation quark q 1 and 
/3e q is the branching ratio of the LQ to I and a quark. For LQ masses exceeding 
y/s, limits are set on the four-fermion contact-interaction term \ eqa A^„/Mlq for 
leptoquarks that couple to an electron and a quark q a and also to £ and a quark 
qp. Some of the limits are also applicable to lepton-flavor-violating processes me- 
diated by squarks in i?-parity-violating supersymmetric models. In some cases 
involving heavy quarks and especially for I = r, the ZEUS limits are the most 
stringent published to date. 
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1 Introduction 



In the Standard Model (SM), lepton flavor is conserved. While the reported observation 
of neutrino oscillations [1,2] implies that lepton-flavor violation (LFV) does occur, min- 
imal extensions to the SM [3] that allow for finite neutrino masses and thereby account 
for neutrino oscillations do not predict detectable rates of LFV at current collider exper- 
iments. However, many extensions of the SM, such as grand unified theories [4], models 
based on supersymmetry [5], compositeness [6] or technicolor [7] involve LFV interactions 
at fundamental levels. 

In high-energy positron-proton collisions at HERA, reactions of the type eg« — > £qj, where 
qi and qf denote initial- and final-state quarks and I denotes a fi or a r with high transverse 
momentum, can be detected with high efficiency and small background. Indirect searches 
for such reactions have yielded very strong constraints [8] for cases where q^ and qf are 
light quarks. However, in some cases involving heavy quarks, especially when i = r, the 
sensitivity of HERA extends beyond existing low-energy limits. 

This paper reports on a search for LFV processes in e + p collisions using data collected by 
the ZEUS experiment from 1994 to 1997 with an integrated luminosity, C, of 47.7 pb" 1 . 
Previous searches for LFV at HERA have been reported by ZEUS [9] ( C ~ 4pb~ x ) and 
HI [10] (C rsj 37pb~ 1 ). 

2 Phenomenology 

There are several mechanisms whereby lepton flavor can be violated in ep collisions. This 
paper considers two main possibilities: leptoquarks and R-parity-violating squarks. 

2.1 Leptoquarks 

Leptoquarks (LQs) are bosons that carry both lepton (L) and baryon (B) numbers and 
have lepton-quark Yukawa couplings. Such bosons arise naturally in unified theories 
that arrange quarks and leptons in common multiplets. A LQ that couples to leptons 
of two different generations would induce LFV. The Buchmuller-Riickl-Wyler (BRW) 
model [11], which assumes the most general Lagrangian with SU(3)c x SU(2) L x U(l)y 
invariant couplings of a LQ to a lepton and a quark, is used to classify LQ species and to 
calculate cross sections for LQ-mediated processes. The following additional assumptions 
were made to simplify the models under consideration: 

1. one LQ species dominates the cross section of the process; 
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2. members of each SU(2) multiplet are degenerate in mass; 

3. LQs couple to either left-handed or right-handed leptons, but not both. 

There are 10 different LQ states in the BRW model, four of which can couple to both left- 
and right-handed leptons. Because of the third assumption above, models in which these 
states have left- or right-handed couplings will be treated separately in this analysis. Each 
state is characterized by spin J = or 1, weak isospin T = 0, 1/2 or 1 and fermion number 
F = or ±2 (where F = 3B + L). Following the Aachen notation [12], scalar (J = 0) 
and vector (J = 1) LQs are denoted Sj. and V^, respectively, where x — L, R denotes 
the chirality of the lepton that couples to the LQ. When two different hypercharge states 
are allowed, one is distinguished by a tilde. In this paper, LQs with couplings \ eqa to an 
electron and a quark q a , and \i qp to a lepton I (/i or r) and a quark qp, are considered^. 
The subscripts a and (3 label the quark generations. The LQ species determines whether 
q a or qp are up- or down-type quarks. In addition to mediating LFV interactions, such 
LQs would also mediate flavor- conserving interactions with an e or a v e in the final state. 
These final states were not searched for in this analysis, but they were taken into account 
in calculating branching ratios (the same is true for final states with or v T ). 

In ep collisions, if the LQ mass, M L q, is below y/s (low- mass LQs), the LQ is predom- 
inantly produced as an s-channel resonance, as shown in Fig. |T|(a). In this case, only 
incident u or d quarks, denoted q±, which couple to the incident positron to produce 
F = LQs, are considered. In the e + p data analyzed here, the production cross section 
for F = LQs is much larger than for F = —2 LQs, assuming that Mlq is sufficiently 
large so that the production is valence-quark dominated, since F = — 2 LQs would be 
produced via e + q fusion. 

For small values of the Yukawa coupling, X eqi , the resonance width becomes negligible and 
the s-channel Breit-Wigner line shape can be approximated (neglecting radiative effects) 
by a 5-function at Mlq = \^xs, where x is the Bjorken variable in deep inelastic scattering 
(DIS). This leads to the Narrow Width Approximation (NWA) 

a™(s, M LQ ) = ( J + 1) ^CV 3 qi (x = Qjj) , (1) 

where T3 is the third component of the weak isospin, Ct 3 is the square of the relevant 
SU(2) Clebsch-Gordan coefficient, and qi(x, Qq) is the valence-quark density in the proton 
evaluated at the scale Ql = M^q. The total production cross section for a given LQ is 
given by the sum over all states of the SU(2) multiplet that couple to a positron and a 
quark. The NWA becomes inaccurate if qi(x, Qq) varies significantly with x on a scale 

1 Note that in the BRW model, some \ — L LQs also have neutrino-quark couplings; these couplings 
are fixed by SU(2)l invariance to be equal to the corresponding charged lepton-quark couplings. 
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corresponding to the LQ width, T L q oc A 2 M L q. In the present analysis, this occurs 
only when Mlq is close to y/s (x — > 1). In this region, q\ falls steeply with x and the 
convolution of q\ with the Breit-Wigner line shape results in contributions to the cross 
section from quarks with x below the resonant peak. These non-resonant contributions 
to the cross section are neglected in the NWA. For LQs that couple to u (d) quarks, 
with M LQ = 270 GeV (250 GeV) and \ eq = 0.3, a NWA underestimates the cross section 
by ~ 20%. The rate for LFV events is proportional to <jflg q where f3z q is the branching 
ratio to the iq final state. In this paper, the NWA is used to calculate cross sections for 
Mlq < y/s, so that limits on a/3i q can be simply converted to limits on X eqi \f]h q . The 
NWA underestimates the cross section, leading to conservative limits. 

If Mlq > y/s (high-mass LQs), both s- and u-channel diagrams contribute, see Fig. [l|. If 
M^q S> y/s, the LQ propagator contracts to a four-fermion contact interaction and the 
cross section is proportional to [A ega A^„ / Mf^] 2 . In this high-mass approximation (HMA), 
the cross section for an F = LQ, in e + p collisions, can be written as 



a 



.HMA 
F=0 



327T 



dxdyxq a (x,s) f(y) + / dxdyxq p {x,-u) g{y) 



(2) 



with 



f{y) 




scalar LQ 
y) 2 vector LQ 



g{y) 



y) 2 /2 scalar LQ 
vector LQ 



where y is the inelasticity, s = sx and u = sx(y — 1) are the scales at which the quark 
densities q a and qp are evaluated. The first and second integrals in (|2|) are due to the s— 
and u — channel contributions, respectively (\F\ = 2 LQs couple a quark in the w-channel 
and an anti-quark in the s— channel). The accuracy of the HMA increases with increasing 
LQ mass. For Mlq > 600 GeV, the minimum mass considered for this high-mass analysis, 
the accuracy is better than 10%. 

In the high-mass case, LQ scenarios are characterized by the 14 LQ species, the three 
generations of q a and q@, and the two possible final-state leptons, leading to a total of 252 
different LQ scenarios. 

NLO QCD corrections [13, 14] were applied only to the NWA production cross section for 
scalar LQs, since no calculation is available for vector LQs or for high-mass scalar LQs. 
These corrections increase the production cross section by ~ 15% at Mlq = 150 GeV, 
increasing to ~ 30% at M LQ = 250 GeV. 

Corrections for QED initial-state radiation (ISR), evaluated using the Weizsacker- Williams 
approximation [15,16], were applied to both the low- and high-mass cases. The QED ISR 
correction reduces the NWA cross section by ~ 3% at Mlq = 150 GeV and by ~ 25% 
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when Mlq approaches the kinematic limit. For high-mass LQs, QED ISR corrections, 
evaluated at Mlq = 600 GeV, were applied. They lower the cross section by less than 5%; 
the corrections decrease at higher masses. 



2.2 H-parity- violating squarks 



Supersymmetry (SUSY), which links bosons and fermions, is a promising extension to 
the SM. It assumes a supersymmetric partner for each SM particle, a bosonic partner for 
a fermion and vice-versa. i?-parity is a multiplicative quantum number defined as R p = 
(_1)3B+l+2J_ For g M part i c i eS; R p = \- for SUSY particles (sparticles), R p = -1. In R p - 
conserving processes, sparticles are pair produced and the lightest supersymmetric particle 
(LSP) is stable. In models with R p violation (R p ), single SUSY-particle production is 
possible and the LSP decays into SM particles. Of special interest for HERA are Rp 
Yukawa couplings that couple a squark (SUSY partner of a quark) to a lepton and a 
quark, which are described in the superpotential by the term [17] X\^ k L l Q' J D , where i, j 
and k are generation indices, L and Q denote the left-handed lepton and quark-doublet 
superfields, respectively and D denotes the right-handed quark-singlet chiral superfield. 
Expansion of the superfields using four- component Dirac notation yields 



C = X' 



ijk 



-e L u L a R 



e L u L Ct 'R 



U, 



(4)* + Wr + v l L d{d R + (j? L y<f L (d 



+ h. c. 

(3) 



The superscript c denotes charge conjugation and the asterisk denotes complex conju- 
gation of scalar fields. For i = 1, the second and third terms will result in and d k 
production in ep collisions. Identical terms appear in the Lagrangians for the scalar lep- 
toquarks Sh 2 and Sq [18]. The coupling gives rise to the reaction e + d — > u 3 L , while 
the coupling X' llk would cause the reaction e~u — > d k R . 

Lepton-flavor violation would occur in models with two non-zero Yukawa couplings involv- 
ing different lepton generations. For example, non-zero values of Ay X and A^- fc (i = 2, 3) 
would yield the process e + d — > v? £ + d k , where i — 2, 3 corresponds to I — //, r. 
Squarks also undergo i? p -conserving decays to a quark and a gaugino, which were not 
considered in this analysis. Low-mass coupling limits on Sy 2 LQs can be interpreted, 
using X eQi \[fitq = X'uiy/ /3&i ->fq, as limits on v? squarks that couple to eq 1 and to £q. In 
the low-mass case, the limits apply for any final-state quark q (except top). High-mass 
LQ limits can also be applied to squarks as described in Section |9.3| . 
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3 The ZEUS detector 



A detailed description of the ZEUS detector can be found elsewhere [19]. A brief outline of 
the components which are most relevant for this analysis is given below. Charged particles 
are tracked in the central tracking detector (CTD) [20], which operates in a magnetic field 
of 1.43 T provided by a thin superconducting coil. The CTD consists of 72 cylindrical drift 
chamber layers, organized in 9 superlayers covering the polar angle^ region 15° < 9 < 164°. 
The transverse-momentum resolution for full-length tracks is (t(j>t)/pt — 0.0058pT © 
0.0065 © 0.0014/pr, with in GeV. The CTD was used to reconstruct tracks of isolated 
muons and charged r-decay products. It was also used to determine the interaction vertex 
with a typical resolution of 4 mm (1 mm) along (transverse to) the beam direction. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. The calorime- 
ters are subdivided into towers each of which subtends a solid angle from 0.006 to 0.04 
steradians. Each tower is longitudinally segmented into an electromagnetic (EMC) section 
and two hadronic (HAC) sections (one in RCAL). Each HAC section consists of a single 
cell, while the EMC section of each tower is further subdivided transversely into four cells 
(two in RCAL). The CAL energy resolutions, as measured under test-beam conditions, 
are a(E)/E = 0.18/VE for electrons and a(E)/E = 0.35/VE for hadrons (E in GeV). 
The arrival time of CAL energy deposits is measured with sub-nanosecond resolution for 
energy deposits above 4.5 GeV, allowing the rejection of non-ep background. 

The FMUON detector [19] consists of layers of limited streamer tubes and drift-chamber 
planes located up to 10 m from the interaction point. The toroidal magnetic fields of 
the iron yoke (1.4 T) that surrounds the CAL and of two toroids (1.6 T) located about 
9 m from the interaction point enable muon momentum measurements to be made. The 
FMUON tags high-momentum muons (muons with momenta below 5 GeV are unlikely 
to emerge from the FCAL) with polar angles in the range 8° < 9 < 20°, extending well 
beyond the CTD acceptance. 

The luminosity was measured by the luminosity detector (LUMI) from the rate of the 
Bethe- Heftier process e + p — > e +r yp [22], where the photon is detected in a lead-scintillator 
calorimeter located at Z = —107 m in the HERA tunnel. The uncertainty on the lumi- 
nosity measurement was 1.6%. 



2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam 
direction. 
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4 Monte Carlo simulation 



The simulation of the LQ signal, including both s- and w-channel processes, was performed 
using the generators LQMGEN 1.0 [23] (low-mass LQs) and LQGENEP 1.0 [24] (high- 
mass LQs) based on the BRW model [11]. Both generators are interfaced to JETSET 
7.4 [25] to simulate hadronization and particle decays. 

The following SM backgrounds were considered: charged current (CC) and neutral current 
(NC) deep inelastic scattering (DIS) were simulated using DJANG06 2.4 [26], with the 
color-dipole model ARIADNE 4.08 [27] used to simulate the hadronic final state. Elastic 
and inelastic 77 £ + £~ reactions were simulated with LPAIR [28]. EPVEC 1.0 [29] 
was used to simulate W production. Photoproduction processes were simulated with 
HERWIG 5.8 [30]. The ZEUS detector and trigger were simulated with a program based 
on GEANT 3.13 [31]. The simulated events were processed by the same reconstruction 
programs as the data. 



5 Kinematic quantities 

Global calorimeter sums were calculated as follows: each calorimeter cell i with an 
energy deposit E^ above a threshold was assigned a four-momentum P M , defined as 
P M = (Ei, EiCos<f)iSm6i, EiSiiKpisinOi, EiCos8i), where (pi and 8i are the azimuthal and 
polar angles of the cell center relative to the event vertex. The total four-momentum 
deposited in the calorimeter (E,Px,Py,Pz) is given by the sum of the four- momenta 
for all cells. The transverse energy, E t , is given by 'Y^ i E i sai9i. The missing transverse 
momentum, P t , is given by a/ P\ + Py- The azimuth assigned to P t , m iss, was defined 
by cos0 m i SS = —Px/Pt and sin0 miss = —Py/P t . Jets used in identifying hadronic r de- 
cays were reconstructed using an (77, 0) cone algorithm [32] with cone radius R — 1. The 
inputs to the jet algorithm were the four-momentum vectors of each calorimeter cell. The 
invariant mass of a jet, Mj et , was calculated from the sum of all four-momentum vectors 
assigned to the jet. The transverse energy of a jet was denoted by E 3 t et . 

The E — Pz of the initial state is twice the positron beam energy, 2E e = 55 GeV. For 
events that are fully contained in the calorimeter (ignoring particles escaping through the 
forward beam hole, which carry negligible E — Pz), the measured E — Pz should be near 
2E e . In photoproduction processes, where the final-state positron escapes through the 
rear beam hole, the E — Pz spectrum falls steeply, so that a cut on E — Pz is useful in 
reducing such backgrounds. Events with high-energy muons, which deposit only a small 
fraction of their energy in the calorimeter, will also have E — Pz substantially below 2E e . 



In the search for the e — » ji transition (see Sections |6.1| and |6.2j), a cut was made on the 
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quantity E — P z + A M , where A^ = f t {l — cos 6^)/ sin^ is an estimate of the E — P z 
carried by the muon, assuming that the transverse momentum of the muon is equal to 
ft, Op is the polar angle of the muon track. 



6 Event selection for the e — ► /x transition 

Events from the reaction ep — > fiX, mediated by a heavy LQ, would be characterized by 
a high-transverse-momentum (P t ) muon balanced by a jet. Since only a small fraction of 
the muon energy is deposited in the calorimeter, these events would have a large f. The 
offline event selection consisted of two steps: a pre-selection of events with f t and a final 
selection requiring an identified muon. 



6.1 Pre-selection 

The trigger, which is identical to that used in the CC DIS measurement described else- 
where [33] , was based on a cut on f with a considerably lower threshold than the selection 
cuts described below. After applying cuts to reject non-ep backgrounds (mainly cosmic 
rays and beam-gas interactions), the following pre-selection requirements were imposed: 

• a reconstructed vertex with Z coordinate |Zvtx| < 50 cm; 

• f t > 15 GeV and f t /^E~ t > 2.5v / GeV; 

• 25 GeV < E - P z + A M < 100 GeV; 

• no electron with energy larger than 10 GeV. 

The third cut discriminates against photoproduction events, while the fourth cut sup- 
presses NC DIS. The electron finder [34] is based on a neural-network algorithm. After 
the pre-selection, 164 events remained, compared with 177.3 ±3.8 events predicted by the 
SM simulation normalized to the integrated luminosity of the data. The error associated 
with the prediction arises from the generated MC statistics. The SM expectation is dom- 
inated by CC DIS, with small contributions from ep — > efi + fi~X and from W production. 
Simulated distributions of E — P z + A^, f t , and f t /\fE~ t are compared with the data in 
Fig. H]. Good agreement is seen. 



6.2 Muon identification 

Two methods of muon identification were employed. The first, for very forward muons 
(8° < 9^ < 20°), required a reconstructed track in the FMUON detector with azimuth 
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within 20° of m i SS . In the second selection, for central muons (15° < 6> /t < 164°), the 
following CAL- and CTD-based requirements were imposed: 

• a track that points to the vertex with transverse momentum (-P 4 trk ) above 5 GeV and 
an azimuth that differs from m i SS by less than 20°; 

• no additional tracks with azimuth within 50° of </> m i SS and P t trk > 1 GeV; 

• the calorimeter energy deposits are consistent with those expected from a minimum 
ionizing particle in an (rj, 0) cone of radius R = 0.3, centered on the track; 

• muons with 115° < 9^ < 130° were excluded to eliminate the background from a very 
small fraction of electrons for which a large fraction of the energy was absorbed in the 
dead material between the BCAL and the RCAL. 

After the muon identification, 2 events are left in the data, while the SM expectation is 
1.43 ± 0.38, mainly from ep — > efi + fi~X. 

6.3 Final selection 

The final selection was designed to reduce the SM background to a very low level. The 
following cuts were applied: 

• P t > 20 GeV; 

• f t /V^>4x/GeV; 

• E - P z + A M > 30 GeV. 

No event survived these cuts, while 0.40 ± 0.18 events are predicted by SM processes, 
mostly from ep — > efi + fi~X. 

7 Event selection for the e — ► r transition 

This channel is characterized by an isolated r with high P t balanced by a jet. Separate 
selections were made for hadronic r decays (65%) and for the leptonic decays r — > iviy T 
(35%). The same trigger as described in Section |6.1| was used. The offline event selection 
consisted of a pre-selection common to hadronic and leptonic r decays and final selections 
specific to each decay mode of the r. These mode-specific selections make use of the fact 
that one or more neutrinos are emitted in r decay producing P t approximately aligned 
with the r. To produce a reasonably large event sample to compare with SM predictions, 
the selections for each r decay mode were done in two steps. 
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7.1 Pre-selection 



In addition to cuts to reject non-ep background, the pre-selection requirements were: 

• a reconstructed vertex with Z-coordinate |2vtx| < 50 cm; 

• 20 GeV < E - P z < 52 GeV; 

• energy in RCAL < 7 GeV. 

The second cut reduces the photoproduction background. The third cut rejects NC DIS 
events where the positron was scattered into the RCAL. 

7.2 Selection of hadronic r decays 

Events with a narrow 'pencil-like' jet consistent with hadronic r decay were selected with 
the following requirements: 

• the transverse energy of the jet associated with the r should satisfy El~ 3e > 10 GeV; 

• Mjet <7GeV; 

• 1, 2 or 3 tracks associated with the jet; 

• the number of calorimeter cells associated with the jet, iV ce ii s , is at least 10 (to suppress 
electrons) and at most 50 (to ensure that the jet is narrow); 

• -^90% < 0.3, where Rgo% is the radius of the (rj, 0)-cone centered on the jet axis that 
contains 90% of the jet energy; 

• /emc < 0.95, where /emc is the fraction of the jet energy deposited in the electromag- 
netic section of the calorimeter; 

• /emc + /lt < 1-6, where /lt is the momentum of the most-energetic track in the jet 
divided by the jet energy (leading-track fraction). 

The last two cuts reject electrons, for which /emc ~ 1 an d /lt ~ 1- After these cuts, 367 
data events were selected in comparison to 377.7 ± 12.5 from the SM expectation (mainly 
from NC DIS, CC DIS and photoproduction). Figure |3| shows several distributions of 
characteristic variables of the r candidates at this stage of the analysis. The SM simula- 
tion provides a reasonable description of the data. 

The final stage of the hadronic r-decay selection requires events consistent with a two- 
body r+jet final state with an invariant mass above ~ 100 GeV: 

• ET iet > 20 GeV; 

• the azimuthal angle of the r candidate is within 20° of m i SS ; 
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• f > 12 GeV; 

• at least one additional jet with E 3 t et > 25 GeV. 

No candidate satisfying these requirements was found, while 0.62±0.18 events are expected 
from SM processes. 



7.3 Selection of r — > ^LV il u T decays 



After the pre-selection described in Section |7.1| , events with an isolated high-P t muon 
balanced by a jet were selected. Isolated muon candidates were identified using a neural- 
network algorithm that analyzed the pattern of longitudinal and transverse energy depo- 
sition in the calorimeter and matching track(s) in the CTD and/or the muon chambers. 
Since the energy deposited in the CAL by the muon is typically a small fraction of the 
energy of the r, cuts on f were applied. The initial requirements were: 

• a muon with P t > 10 GeV and 8° < 0„ < 125°; 

• f > 15 GeV; 



• ft/VE~ t > 4 v / GeV; 

• a jet with E\ ct > 25 GeV; 

• events with an identified electron [35] with energy greater than 10 GeV were vetoed. 

The last cut reduces NC DIS background. Figures §(a-b) show the distributions of E — Pz 
and f for the r — > [ivv candidates after these cuts, compared to the SM background. 
Good agreement is observed. After these cuts, 119 data events remained, compared to 
107.2 ± 7.4 events from the SM expectation (mainly CC DIS and photoproduction). 

The final selection consisted of two cuts: 

• f > 20 GeV; 

• the muon azimuth differs from m i SS by less than 20°. 

No event passed the final selection, while 0.23 ± 0.07 events were expected from SM 
processes. 



7.4 Selection of r — > eu e u T decays 

After the pre-selection described in Section |7.1| , events with an isolated electron and a jet 
were selected by imposing the following requirements: 

• an identified electron [35] with energy greater than 20 GeV and the polar angle 9 e 
satisfying 8° < 6 e < 125°; 
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• f t > 10 GeV; 

• a jet with Ef * > 25 GeV. 

Figures f§(c-d) show the distributions oi E — P z and f t for the r — > eVu candidates that 
satisfied these requirements, where 116 data events were selected and 109.1 ± 5.4 events 
from SM backgrounds were expected (mainly NC DIS and photoproduction). 

The final selection consisted of a higher f cut and a requirement that the lepton and the 
jet be back-to-back: 

• ft > 15 GeV; 

• the azimuth of the electron differs from m i SS by less than 20°. 

No event passed the final selection, while 0.32 ± 0.10 events were expected from SM 
processes. 



8 Efficiencies 

The selection efficiencies were evaluated using signal MC events (see Section For res- 
onant production of lepton-flavor-violating scalar LQs, the /i-channel selection efficiency 
falls from 60% to 52% as Mlq increases from 140 GeV to 280 GeV, while the efficiency for 
vector LQs drops from 64% to 56%. For M LQ > 240 GeV, the FMUON-based muon se- 
lection increases the selection efficiency by about 20% compared to the CAL-CTD-based 
selection alone. Over the Mlq interval from 140 GeV to 280 GeV, the selection efficiency 
for LQs that couple to r increases from 24% to 31% for scalar LQs and from 21% to 33% 
for vector LQs. 

For LQs with M L q ^> y/s, the efficiencies are almost independent of Mlq, but depend 
strongly on the generation of the initial-state quark. For e — > \x transitions, the selection 
efficiency ranges from 15% to 45% for F = LQs and from 15% to 35% for \F\ = 2 LQs. 
For e — ► r transitions, the efficiencies are lower and range from 5% to 19% for F = 
LQs and from 4% to 16% for \F\ — 2 LQs. When the initial-state quark is a sea quark 
and especially for s, c, or b quarks, the efficiency is considerably lower than for valence 
quarks due to the softer x spectrum, which results in a lower transverse momentum for 
the final-state lepton. 
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9 Results 



Since no candidate for LFV processes was found, limits were set on these processes. All 
limits were evaluated at 95% CL. using a Bayesian approach, assuming a flat prior for 
the signal cross section. Systematic uncertainties in the detector simulation and in the 
integrated luminosity (see Section |9.1|) were taken into account using a method described 
elsewhere [36]. For low-mass LQs with narrow width, the branching ratio, (3, was regarded 
as a free parameter and limits were set on a Peg- These limits were converted to limits on 
K qi \ jlhq using (jip corrected for QED-ISR and NLO QCD (only for scalar LQs). For high- 
mass leptoquarks, the cross-section limit was converted to a limit on X eQa \i qf3 / M1q using 
(|2|) with QED-ISR corrections. The CTEQ4 [37] parameterizations of parton densities 
were used to evaluate cross sections. 



9.1 Systematic uncertainties 

The uncertainty on the integrated luminosity is 1.6%. Systematic uncertainties of 3% on 
the CAL energy scale and 10% on the CAL response to muons were taken into account. 
The resulting variations on the efficiency for the muon (tau) channel were 3% (4%) for 
low-mass LQs and up to 15% (17%) for high-mass LQs that couple to b quarks in the 
initial state. 

Systematic uncertainties in the cross-section evaluation, related to the choice of parton 
density function (PDF), were investigated using MRST [38] as an alternative choice to 
CTEQ4. The main differences were found for low-mass LQs with masses close to y/s when 
very high-x quarks are involved. In these cases, limits calculated using MRST were stricter 
than the CTEQ4-based limits presented here. Another possible source of uncertainties 
for vector and high-mass scalar LQs are the unknown NLO-QCD cross-section corrections 



(see Section 2.1 



9.2 Low-mass LQ and squark limits 

Figure ^| shows upper limits on oft m and a[3 Tq . For e — > fi, the search is sensitive to 
processes with cross sections as low as 0. 1 pb, while for e — > r, the sensitivity is 0.2 - 0.3 pb. 
These limits apply generally to narrow resonances with LFV decay modes, for example, 



to the tfl p squarks described in Section [272. 

Upper limits on \ eq \ffh q have been derived for F=0 LQs by assuming resonantly pro- 
duced LQs described by the BRW model. These limits can be applied to processes in- 
volving any quark generations in the final state (excluding the t— quark). Figures |6](a-b) 



12 



show the upper limits on X eq ■sj & m for scalar and vector LQs, respectively. Under the 
assumption that \ eq ^ = X M0 , limits on X eqi can be derived. These are compared to limits 
from low-energy experiments in Figs. |6](c-d) for Sy 2 and Vq 1 LQs. These states do not 
couple to neutrinos and therefore (3 m = 0.5. For Mlq < 250 GeV, the ZEUS limits are 
stronger than the low-energy limits for LQs that couple to /i and b. Limits on X eQi (3 Tq 
are shown in Figs. |7|(a-b). Figures 0(c-d) show the corresponding limits on X eq , assuming 
A e?i = X Tqf3 (and therefore (3 Tq = 0.5). The ZEUS limits are more stringent than the limits 
from low-energy experiments over a wide mass range, with the exception of limits from 
K + — > tx+vv [39]. As described in Section |2.2[ , the limits on X e q 1 ^/Peq f° r Sy 2 can be 
interpreted as limits on \\ji^J flvi^iq for v? squarks. 

Another way to illustrate the sensitivity is to assume that the couplings have electromag- 
netic strength (X eg = Xi qp = 0.3 ~ y/Aira). In this case, LQs with masses up to 283 GeV 
are excluded, as shown in Table [l]. Alternatively, as shown in Table 0, for a fixed Mlq of 
250 GeV, values of X eqi ^J~fk q down to 0.020 (0.027) for LQ -> /iq (LQ — > rq) are excluded. 

The CDF [40] and D0 [41] collaborations exclude scalar LQs with M LQ < 202 GeV 
and M LQ < 200 GeV, respectively, at 95% C.L. with $ m = 100%. CDF [42] excludes 
-Mlq < 99 GeV with j3 Tb = 100%. The ZEUS limits are complementary to those of the 
Tevatron in the sense that the latter are independent of the Yukawa couplings and assume 
that LQs couple to a single lepton generation. 

9.3 High-mass LQ and squark limits 

For M L q ^> y/s, limits on \ eqa \i qfj / were evaluated for all combinations of quark 
generations (a, (3). Tables |3| and |] show these limits for F = and \F\ = 2 LQs, 
respectively, that couple to fxqp. Tables [5] and |6| show the corresponding limits for the 
LQs coupling to rqp. In many cases involving c and b quarks, the ZEUS limits improve on 
the low-energy limits [8,43,44]. Limits obtained by HI [10] are comparable to the ZEUS 
limits. 

The limits on A^A^/Mlq for Sh 2 can be interpreted as limits on A^^A-^/M? for a 
w-type squark of generation j, where i = \x or r for % = 2 or 3, respectively. Similarly, 
the limits on Sq LQs can be interpreted as limits on A' lafc A^ A ./M| for a d-type squark of 
generation k. 

10 Conclusions 

A search for lepton-flavor violation has been performed with 47.7 pb -1 of e + p data at 
v/s = 300 GeV collected with the ZEUS detector at HERA in 1994-1997. Both the /i and 
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r channels have been analyzed. No evidence for LFV processes has been found. 

Limits at 95% CL. on cross sections, couplings and masses for F=0 LQs that mediate 
LFV processes have been set. Assuming the couplings X eQi = \e qp = 0.3, lower mass limits 
between 258 and 283 GeV have been derived for various LQs decaying to fiq or rq. For 
M L q = 250 GeV, upper limits for \ eq \J fi m in the range (2.0 - 10) • 10~ 2 and for A e?i ' f3 Tq 
in the range (2.7-15) ■ 10~ 2 were obtained. Limits on Sy 2 also apply to up-type squarks 
that have i?-parity-violating couplings to both a positron and either a \i or a r. 

For LQs with Mlq \fs, upper limits on X eqa Xe qi3 / M 2 q have been obtained and compared 
with bounds from low-energy experiments. Some of these limits also apply to high-mass 
R p - violating squarks. A number of ZEUS limits are the most stringent published to date, 
especially for e — > r transitions. 
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LQ type 


oL 
°l/2 


CL 
°l/2 


CR 
°l/2 




K R 


% R 


V? 


//-channel limit on Mlq [ GeV] 


263 


278 


278 


261 


266 


280 


283 


r-channel limit on Mlq [ GeV] 


258 


275 


276 


259 


263 


277 


282 



Table 1: The 95% CL. lower limits on M L q for the /i- and the r-channels 
assuming \ eqi = \i q = 0.3. 



LQ type 


CL 
°l/2 


CL 
°l/2 


cR 
°l/2 


v L /v R 


V* 




/i-channel limit on \ eqi fi m 


0.10 


0.038 


0.036 


0.081 


0.029 


0.020 


r-channel limit on A e(3i ' (3 Tq 


0.15 


0.054 


0.051 


0.10 


0.038 
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Table 2: The 95% CL. upper limits on A e9i y^ /or a leptoquark with mass 
Mr n = 250 GeV. 
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Table 3: Upper limits at 95% C.L. on Ae^A^/M]^ in units ofTeV~ 2 , for F = 
LQs that couple to eq a and to \iqp. The columns correspond to the F = LQ 
species. The eq a combination for the s-channel case is reported under the LQ type. 
Each row corresponds to a different combination of quark generations (a, (3) which 
couple to the positron and the \x, respectively. Within each cell, the measurement 
which provides the most stringent low-energy constraint is shown on the first line 
and the corresponding limit [8,43,44] ^ s given on the second line. The ZEUS limits 
are shown on the third line of each cell ( enclosed in a box when stronger than the 
low-energy constraint). The * indicates cases where a top quark must be involved. 
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Table 4: Upper limits at 95% C.L. on Ae^A^/M^g m units o/ TeV~ 2 , for 
\F\ = 2 LQs that couple to eq a and to fiqp. The columns correspond to the \F\ = 2 
LQ species. The format of the table is described in the caption of Table^. 
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Table 5: Upper limits at 95% C.L. on \ eqa \ Tqf J M"Iq in units of TeV~ 2 , for 
F = LQs that couple to eq a and to rqp. The columns correspond to the F = 
LQ species. The format of the table is described in the caption of Table f|. 
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Table 6: Upper limits at 95% C.L. on \ eqa \ rqf J M"Iq in units of TeV -2 , for 
\F\ =2 LQs that couple to eq a and to rqp. The columns correspond to the \F\ =2 
LQ species. The format of the table is described in the caption of TableQ 
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Figure 1: (a) s-channel and (b) u-channel diagrams contributing to LFV processes 
induced by F = LQs. In e + p scattering, \F\ = 2 LQs couple to antiquarks in the 
s-channel and to quarks in the u-channel. 
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Figure 2: Distributions of event variables after the ^-channel pre- selection for 
data (solid points) and SM simulation (shaded histograms) for a) E — Pz + A M , b) 
P t and c) f t /y/E~ t . The dashed histograms simulate the signal from a scalar LQ 
with Mlq = 260 GeV normalized to the 95% C.L. cross-section upper limit (see 
Section 
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Figure 3: Comparison of data (solid points) and simulated SM background (shaded 
histograms) for candidate jets from hadronic r decays. The distributions are dis- 
played for events that pass the selection cuts described in Section \7.3j except the 
ones imposed on the variable considered (indicated by the arrows). Shown are the 
distributions of: (a) 9 T .j et , the polar angle of the r candidate jet; (b) E T t ~ 3et , the 
transverse energy; (c) iVceiis, the number of calorimeter cells belonging to the jet; 
(d) Rgo%> the (r], <p) -radius containing 90% of the jet energy; (e)fEMC, the fraction 
of the jet energy in the EMC section of the calorimeter; (f) /emc + /lt where /lt is 
the momentum of the leading track divided by the jet energy. The SM backgrounds 
include NC and CC DIS, photoproduction, W production and 77 — > t + t~ . The 
dashed histograms simulate the signal from a scalar LQ with a mass of 260 GeV 
normalized to the 95% C.L. upper limit on the cross section (see Section^). 
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Figure 4: Comparison of data (solid dots) with simulated SM background (shaded 
histogram) for the distributions of a) E — Pz and b) P t for the t — > fiuu selec- 
tion. The same distributions for the t — > evv selection are shown in c) and d), 
respectively. The SM backgrounds include NC DIS, photoproduction, CC DIS, W 
production and 77 — > t + t~. The dashed histograms simulate the signal from a 
scalar LQ with a mass of 260 GeV normalized to the 95% C.L. upper limit on the 
cross section (see Section^). 
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Figure 5: The 95% C.L. upper limits on a/3i q as a function of Mlq for scalar 
(full line) and vector (dashed line) LQs for a) LQ — > fiq and b) LQ — > rq. 
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Figure 6: Upper limits on \ eq ft m vs. Mlq for a) scalar and b) vector LQs. The 
quark flavors that couple to the LQs in the initial state are shown in parentheses 
following the LQ species. Upper limits on \ eQi under the assumption (3 m = 0.5 
are shown in c) for scalar LQs and d) for vector LQs that couple to d-type quarks. 
Also shown are existing limits [8,43,44] (dashed lines). The numbers in parentheses 
indicate the generations of the quarks that couple to the e and the \x, respectively. 
The regions above the curves are excluded at the 95% C.L. 
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Figure 7: Upper limits on X £ 



vs. 



^lq for a) scalar and b ) vector LQs. The 



quark flavors that couple to the LQs in the initial state are shown in parentheses 
following the LQ species. Upper limits on X eqi under the assumption (3 Tq = 0.5 
are shown in c) for scalar LQs and d) for vector LQs that couple to d-type quarks. 
Also shown are existing limits [8,43,44] (dashed lines). The numbers in parentheses 
indicate the generations of the quarks that couple to the e and the r, respectively. 
The regions above the curves are excluded at the 95% C.L. 
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